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Abstract: It has been argued that the Swiss-Cheese cosmology can mimic Dark Energy, 
(H when it comes to the observed luminosity distance-redshift relation. Besides the fact that 

Q-i this effect tends to disappear on average over random directions, we show in this work that 

based on the Rees-Sciama effect on the cosmic microwave background (CMB), the Swiss- 
Cheese model can be ruled out if all holes have a radius larger than about 35 Mpc. We also 
show that for smaller holes, the CMB is not observably affected, and that the small holes 
can still mimic Dark Energy, albeit in special directions, as opposed to previous conclusions 
in the literature. However, in this limit, the probability of looking in a special direction 
where the luminosity of supernovae is sufficiently supressed becomes very small, at least in 
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the case of a lattice of spherical holes considered in this paper. 
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1. Introduction 

Ever since it was observed that distant supernovae (SN), of type la, appear dimmer than 
expected in a matter dominated spacetime [1, 2], when combined with measurements of 
the local Hubble factor [3] , cosmologists are led to the conclusion that the recent expansion 
of the universe definitely accelerates. Observations of the Cosmic Microwave Background 
(CMB) prefer a large angular diameter distance to z = Zdec ~ 1100, which in ACDM 
indicates a close to spatially flat universe [4, 5, 6, 7, 8]. Given constraints on the matter 
content of the universe, from CMB and Large Scale Structure, this is achieved when a 
Cosmological Constant is present, given the assumption that the locally observed Hubble 
factor is equal to the global Hubble factor. The amount of clustering of galaxies is slowed 
down by accelerated expansion with respect to a pure cold dark matter universe [9]. The 
latest addition to this impressive evidence are the observed Baryon Acoustic Oscillations 
(BAO), which provide an observation of one length scale at different redshifts, and thereby 
measure the expansion over different times [10, 11]. 

Presented as such, the evidence for acceleration is convincing, see Ref. [12] though for 
critique. However, a number of assumptions needs to be made in order to come to the 
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conclusion that a Cosmological Constant or a form of Dark Energy (DE) is there. If the 
locally observed Hubble factor cannot be extrapolated to the global expansion rate, as is 
the case if, e.g., the observer lives in a large void, the observed angular diameter distance 
of SN may be explained without DE [13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 
27, 28, 29, 30, 31, 32, 33, 34]. Spatial flatness would in that case be achieved with a matter 
density equal to the critical density, Qm = 1, fitting the CMB [34]. A globally small Hubble 
constant is needed in that case to explain the age of the universe [35] . The clustering of 
galaxies can also be caused by a hot dark matter component, such as neutrino's with a 
mass of m u ~ 0.5 eV [34]. The BAO data are obtained as interpreted data, in the sense 
that one has to apply a fiducial cosmology to the data in order to abstract the BAO. To 
our knowledge this has only been done assuming an FLRW-universe, and it is unclear how 
the BAO would be affected by applying different fiducial models. 

There is an ongoing debate as to whether today's universe is properly described by 
a perturbed FLRW-metric, assuming large scale homogeneity and isotropy, or not. See 
Refs. [36, 37] for explications why the usual argument, that the universe today is everywhere 
described by a small Newtonian potential, may not be sufficient to rule out the role of 
inhomogeneities. If the universe is properly described by the FLRW-metric, then the 
evidence for DE is compelling. If it is not, then we need to understand why we observe an 
apparent acceleration, which is phenomenologically well described by ACDM. Many works 
have been devoted to possibilities linking the apparent acceleration to structure formation 
and to a departure from homogeneity of the Universe, none of them however convincing 
the community that structure formation explains the observations without the need for 
DE. See Refs. [38, 39, 40, 41, 42] for reviews. 

Here we will focus on one such proposition: the Swiss-Cheese Cosmology [43, 44, 45, 
46, 47, 48, 36, 49, 50, 51, 52, 53, 54, 55, 56, 57]. The Swiss-Cheese Cosmology is described 
globally by an FLRW-metric, but locally contains holes (voids surrounded by a mass shell) 
described by a Lemaitre-Tolman-Bondi metric (LTB). The holes match flawlessly to the 
FLRW-metric at the borders, and are mass compensated in such a way that, from outside 
such a hole, the global effect of a patch containing such a hole is as if it were FLRW 
(cheese), due to Birkhoff's theorem [65]. The Swiss-Cheese 1 toy model is an exact solution 
to the Einstein equations, hence does not suffer from averaging problems. Its goal is to 
approximate today's inhomogeneous universe in an exactly solvable manner, with most 
of the matter in structures, separated by voids. Recently it was shown that if curvature 
is a function of space, large local curvature at low redshift, as in this model, is hardly 
constrained by observations [58]. 

The general idea of the Swiss-Cheese Cosmology is that photons travel through holes 
and structures, where holes have a lensing effect such that distant supernovae, observed 
through a number of holes, appear dimmer than in a homogeneous EdS Universe. This 
effect has been explored in the literature and it was found: 

• that if holes are perfectly aligned, a lensing effect can mimic a DE of about Ode = 0.4, 

1 The author must express his doubts about the name 'Swiss Cheese'. The name implies that only swiss 
cheese contains holes, and at the same time that all swiss cheese contains holes. 
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• however, that averaged over many random distributions of holes, the effect vanishes, 
given that the holes are spherical. 

In this work, we will show that, besides the on average disappearing of the lensing effect, 
the Swiss-Cheese Cosmology leaves a significant imprint on the CMB. It has been claimed 
that the effect that the holes have on the redshift of photons is marginal in comparison with 
the redshift at which the distant supernovae are observed. This reasoning does hold, but the 
effect is significant when compared to the observed anisotropies in the CMB-temperature, 
which are of the order O(10~ 5 ). For the first time we will show full sky CMB-maps 
in a Swiss-Cheese Cosmology, for different realisations, each realization with a different 
constant size of holes, ry. When showing these maps, we neglect the primordial power 
spectrum of the CMB-anisotropies, as the secondary anisotropies in most of the cases are 
overwhelmingly larger. These anisotropies are due to the Rees-Sciama [59, 60, 61] effect, a 
non-linear late integrated Sachs- Wolfe effect due to structure formation. 

At the same time, we will show full sky maps of the angular diameter distance, d\, at 
a fixed redshift in each direction. We confirm the results of Ref. [51], that averaged over all 
directions, the angular diameter distance will appear as if the universe were exactly EdS, at 
least when the holes are spherical. This leads to the conclusion that a Swiss-Cheese model 
that has predominantly holes with a radius larger than 35 Mpc is ruled out, since it either 
cannot explain the supernovae in every direction, or it will leave a significant Rees-Sciama 
imprint on the CMB that is already ruled out by observations. 

Our conclusion only applies to cosmologies in which the universe contains a very large 
number of large voids. In principle a smaller number of large voids is not ruled out by 
our analysis. In such a case, however, the voids would play no role at all in explaining the 
distant supernovae, which is not part of the scope of this paper. 

We find that, in the density profile we consider, the maximal effect that the structure 
can have of the luminosity-distance-redshift relation is only marginally dependent on the 
size of the holes. This is a different conclusion than the one drawn in Refs. [36, 53]. We 
will briefly address this difference. 

The different dependence on the size of holes, for different choices of density profile, 
leaves the door open to an 'Apollonian Gasket-like configuration' 2 , in which the universe 
on average is FLRW, but locally is FLRW nowhere. In that case, an observer would see 
through holes in all directions. The CMB, if the configuration is such that all holes are 
sufficiently small, would be left in agreement with established perturbation theory. The 
holes can not be spherical, though, as the cancellation on average of the lensing effect [51] 
still applies to tiny holes as well. Probing such a configuration goes beyond the scope of 
this work, as the model would need significant changes: the holes must be typically smaller 
than rhi ~ 35 Mpc, they must not be spherical, and they must be such that in all directions 
the chance of looking through holes is higher than looking through cheese. Altogether this 
poses serious difficulties for the Swiss-Cheese model. Besides, one might question if such a 

2 The Apollonian Gasket is a fractal, constructed by starting out with three tangent circles, and subse- 
quently filling areas between those circles with circles tangent to the previously drawn circles. See Ref. [62] , 
and figure 1 in [48]. 
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configuration approaches nature in any way, and if the effect remains significant with more 
natural density profiles. 

In section 2 we will briefly overview the metric and related equations describing the 
Swiss-Cheese cosmology and we will present the techniques used to calculate the full-sky 
maps of temperature and angular diameter distance. In section 3 we present the full-sky 
CMB in a Swiss-Cheese universe, for different typical hole sizes. In section 4 we will discuss 
the overall effect of the Swiss Cheese on distance measures. We conclude in section 5. 



2. The Model 



2.1 The metric and geodesies 

The metric We define the Swiss-Cheese model identical to the model in Ref. [48]. In 
the Swiss-Cheese model, the metric is anywhere of the LTB form 

y ,2 (r, t) 



ds z 



-dt 2 + 



W 2 {r) 



-dr z +Y\r, t)dn 2 



(2.1) 



where Y'(r, t) = d r Y(r,t). We can choose the coordinate system such that the universe is 
divided in equally sized cubic boxes, with sides of length 2r c h, each with the origin of its 
own coordinate system in the center of that box. The time coordinate t is scaled such that 
to = corresponds to today, and t = — 1 corresponds to the time of the Big Bang. 

In the cheese, we recover the FLRW-metric by setting Y c ^(r, t) = ra c ^(t) and W^Jr) = 
1 — kr 2 . In the holes, the metric is determined by 



Y hl (r,t) 



M(r) 



3^1hi(r,t) 



with initial conditions at time t 



-0.8, 



W£(r) - 1 



= E(r) 
M(r) 

p(r) 
Y(r,t) 



1 



-H 



2 

E47T 



FLRW 
r 



+ 2E(r), 



1 M(r) 
6ir r 



p(u)Y hl (u, t)Y^(u, t)du, 



o 

Ae 2 CT 2 



Pch 



+ e 



for r < Thi, 
for r > Thi, 



-r. 



(2.2) 

(2.3) 
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In this setup, the free parameters are those describing the matter distribution in the 
holes: the size of the spherical LTB-metric rhi, the minimum density e, the energy den- 
sity scale A, the comoving radius at which the peak of the mass distribution resides ryi, 
the comoving width of the mass shell a, and the FLRW-energy density p c h- The mass 
distribution is constant in time, and a function of comoving radius r only. In physical 
coordinates, however, mass will be moving outwards to form a shell close to the border of 
the LTB-patch. 
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The LTB-evolution is 'switched on' manually at time t = —0.8, such that at — 1 < t < 
—0.8 spacetime is described by the homogeneous and isotropic FLRW-metric everywhere. 
In this way, the energy density and the scale factor are described by continuous functions of 
time everywhere. The spatial curvature W(r), Eq. (2.3), however, appears instantaneously 
at time t = i. This discontinuity plays no role in the quantities investigated here, henceforth 
it is taken for granted in this toy model. 

For the derivation of and motivation for the choices made above, we refer the reader 
to Ref. [48]. 

Geodesic equations As a consequence of the spherical symmetry of each coordinate 
patch, any geodesic will lay in a spatial plane. Hence, without loss of generality, we can 
write the geodesic equations as four independent equations, as in Refs. [47, 48], 

-4^jy *(0)=0, (2-7) 

t(0) =0, (2.8) 
r(0)=r obs , (2.9) 




Y 2 (r, t) 



<K0) =<pobs, (2.10) 



where r ODS and OD s define the location of the observer. 

The constant is defined by the physical angle between the photon geodesic and a 
radial geodesic pointing (which itself has c$ = 0) towards the photon geodesic. At any 
time in any point, the constant satisfies the relation, 

cosa = g l (y)x i (y)g ij (y), (2.11) 
C(j) = (1 + z)Y(r,t)sma, (2.12) 

where g denotes the spatial direction of the geodesic at point y, and x denotes the spatial 
part of a radial geodesic pointing to coordinate y. 

2.2 Dimensions and configurations 

We chose the parameters as in Ref. [48], being r c h = = 0.042k, e = 0.0025, a = rhi/10, 
ryi = 0.037k, A = 50.59 and pflrw = 25, with however a freedom to chose the rescaling 
factor k. For any k, in physical dimensions these numbers correspond to holes with a radius 
of 350k Mpc, in cubes with sides of length I = 2r c h = 700k Mpc, and 5/k holes between 
the observer and t = —0.8 in an optimal direction, as illustrated in figure 1. The cheese is 
chosen to be spatially flat, with 0^=1, i.e., EdS. We use units in which c = 167tGn = 1. 
The translation from these dimensions to physical dimensions is given in Table 1. Given 
the freedom to rescale, the choice of normalization is arbitrary. Throughout this work 
we will always have one size for all holes on a regular lattice within each realization of a 
Swiss-Cheese universe. Secondly, we will always place the observer in the cheese, at a spot 
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Quantity 



Notation 



Unit 



Value 



mass density 
time 

comoving radial coordinate 
metric quantity 
expansion rate 
spatial curvature term 



p(r,t), p(r,t) 
t, T, t, tBB, To 
r 

Y(r,t) 
H(r,t) 
W(r) 



Pco 

(6irpco)- 1/2 
(6irp C0 )- 1/2 
(6ttpco)- 1/2 
o) 

1 



(67rp c ^ 1/2 



9.2 x 1CT 30 g cm" 3 
9.3 Gyr 
2857 Mpc 
2857 Mpc 

3 
2 



3 Hq Qbs 



Table 1: The units used throughout this work, where c = 167tGn = 1- The present critical density 
is pco = SHq 2 obs /8TT, with H a . bs = 70 km s _1 Mpc -1 . Table taken from Ref. [48] 
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Figure 1: A two-dimensional illustration of the configuration of holes throughout the cheese. Gray 
corresponds to the FLRW-mctric, white corresponds to the LTB-mctric. The difference between 
the left and the right illustration is Ki c ft = 2K r i g ht- Throughout this work we will always have one 
size for all holes within each realization of a Swiss-Cheese universe. 



where two adjacent holes are closest to each other. There is no particular reason to choose 
this location, other than for simplicity, and this choice is not important for the conclusions. 

2.3 Different methods for the angular diameter distance 

The angular diameter distance, d\, is defined as the ratio of the size of an object at some 
distance of an observer, and the angular diameter at which it is observed. When, in the 
configuration of aligned bubbles, an observer is located in the cheese at a point where two 
adjacent bubbles are closest to each other, and the observer sees a source for which the 
geodesic connecting the observer and the center of the source is a straight line exactly 
through the centers of a number of bubbles, the angular diameter distance to that source 
amounts to 

4(A) = nrW.«(A))d.«W (2 . 13) 
a 

with a the small angle between the radial geodesic and a geodesic from an edge of the 
source to the observer, the latter geodesic described by r(A), t(X) and 0(A) [48]. For small 
a this quantity is independent of a. In an actual (numerical) calculation, one can integrate 
the geodesic equations for successively smaller alpha, until the quantity converges to 
the same number along the geodesic. 

Eq. (2.13) obviously only holds for the special case of the aforementioned purely radial 
geodesic connecting source and observer. The same method could be used in any direction, 
shooting two photons in almost but just not the same direction. This method would be 
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computationally expensive though, since for each direction one has to reintegrate several 
times the geodesic equation, at very high accuracy such that the relative error in the 
tiny displacement orthogonal to the direction of the photon is small. The orthogonal 
displacement itself will already be extremely small compared to the distance travelled 
along the geodesic, given the small angle at which one shoots. 

In stead, equations describing the exact beam size along the geodesic are, as in [47], 

dO 2 , . N / dt \ 2 ^ 2 



1 d 2 VA 2 , N ( dt\ 2 9 . 

"jr ~° ■ ( 2 - 16 ) 



y/A d\ 2 3 ry ' ' \dX) 



Here A denotes the beam size, 9 is the beam expansion, defined through 9 = fx > an< ^ a 
denotes the beam shear. The beam stretching becomes 9 ± a in two orthogonal directions. 
We refer the reader to Refs. [47, 63] for further explanations. The angular diameter distance 
e^A and the luminosity distance g?l to an observer at any point along a geodesic pointing 
from / to the observer then become, 



d A {\) = J£$- (2.17) 



'source 
\2 



d L (A) = (l + *rd A (A), (2.18) 

with ^source the solid angle at which the source is observed. The initial conditions for a 
beam are, 

dVA 



dX 
r A 



— \Z^sourcej (^"^^) 

A=0 

= 0. (2.20) 

A=0 

With these initial conditions, 9 and a are ill-defined at the initial conditions, as lim^o = 
00. If we write £ = Aa, we find the only relevant equations, 

2. rK n /,u\ 2 



dWA 2 r- . fdty g 



A p(r,t)(-^ (2.22) 



dX 2 3 ry ' ' \dX ) A 3 / 2 

with the extra initial condition £|^_ = 0, if we demand that a is well behaved, i.e. finite, 
at all times. With these equations, the calculation of the angular diameter distance, and 
thereby the luminosity distance, amounts to simply including two more equations in the 
(numerical) integration scheme. Since d\, = (1 + z) 2 d A , hereafter we will mainly focus on 
d\ for simplicity. 

We checked that both the intuitive, Eq. (2.13), and the exact, Eqs. (2.21, 2.22), method 
agree in the case where both can be applied. 
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2.4 Using 2D geodesic equations in a 3D setup 

The spatial plane in which the photon's geodesic lays, is defined by the spatial part of 
a radial geodesic and the spatial part of the photon's geodesic. In the practice of this 
calculation, a photon exits one coordinate patch (A) and enters another patch (B) when one 
of its cartesian coordinated reaches r c h- Let g denote the spatial direction of the geodesic in 
cartesian coordinates, let Xq denote the cartesian coordinates at which the geodesic entered 
patch A, and let v = v/\v\ for any spatial vector v. In cartesian coordinates, the direction g 
is invariant under translations, hence invariant under a transformation from one coordinate 
patch to another. This simplification only holds because the transformation is done at a 
point where spacetime is described by the spatially flat FLRW-metric. During a time step 
of integration of the geodesic equation, corresponding to photons leaving patch A, whilst 
moving in a plane spanned by the orthonormal vectors ef 1 and e^, the transformation at 
the border at time A cross is done as follows: 

<? B (Across) = 9 A ( A cross ) , (2.23) 

= + ef r A cos A + e£rA sin A ± I , (2.24) 




~B\~B aB 



(2.25) 



r a = \a%\, (2.27) 

(2.28) 

cosa B = 5 V B <7§, (2.29) 
(1 + z)Y(r,t) sin a B . (2.30) 



B 



where all quantities are evualated at A croS s- The integration then continues in box B with 
initial values r B and Bj defined in the plane spanned by the orthonormal vectors e B and 
ef . In this example we assumed that it is the third cartesian coordinate of vector g that 
hits the border of the patch. Figure 2 gives a visual explanation of the relation of different 
vectors. 



3. The CMB 

3.1 Temperature maps and their power spectra 

As the universe is described by the matter-dominated FLRW-metric at — 1 < t < t = —0.8, 
photons obtain no anisoptropies from the Swiss-Cheese structure in that period. It suffices 
to integrate photon geodesies backwards in time from today back to t = t, and then to 
compare redshift a photon has experienced with the redshift experienced on average in all 
directions. This translates directly to a relative temperature difference, as T tx a -1 oc 1+z. 
Hence, shooting photons in all directions and obtaining the redshift at time t = —0.8 as 
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Figure 2: A cartoon of a geodesic in two adjacent patches. The geodesic is always in a plane 
in each patch, spanned by orthonormal vectors e\ and e.%, both defined when entering a patch. 
At the exact border of the coordinate patches we have in three-dimensional cartesian coordinates 

cross ) = 9 (Across)- 




Figure 3: Full-sky maps of secondary CMB-anisotropies induced by the Swiss-Cheese structure, 
in Mollweide-projection. The quantity shown is T ~^ T . These temperature maps would be observed 
in Swiss-Cheese universes with holes of size ry = 3.5 Mpc (upper left figure), fhi = 350 Mpc (upper 
right figure), ru = 875 Mpc (lower left figure) and rm = 1.75 Gpc (lower right figure). The cold 
spots correspond to photons that experienced a large Rees-Sciama effect, hot spots correspond 
to photon's that traveled through a non-integer number of holes, and all ring-like structures are 
artifacts of the regular distribution of holes in this configuration. 

a function of direction, one obtains a temperature anisotropy map of the CMB with only 
secondary anisotropies, caused solely by the Swiss-Cheese matter distribution. 

We show full-sky CMB anisotropy maps for four different hole sizes in figure 3. The 
quantity shown is T ^ T , where T is an average over all directions. In each direction, the 
photon departed at the same time t = i, towards us. The hole sizes for these maps are 
rhi = 3.5 Mpc (top left), ry = 350 Mpc (top right), = 875 Mpc (bottom left) and 
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Figure 4: The C\ • of secondary anisotropics for different Swiss-Cheese cosmologies, in bins of 5 
multipoles. For comparison we plot the full C ; TT -spectrum in a standard ACDM cosmology with 
^baryon = 0.045, Q c dm = 0.245, = 0, ^de = 0.71, h = 0.7 (red solid line, unbinned). A Swiss- 
Cheese universe with holes of radius smaller than 35 Mpc, potentially leaves the CMB unaltered 
with respect to the standard cosmological model. Note the numerical limitations beyond the ankle 
at high I ~ 100, as discussed in the text. 



r^i = 1.75 Gpc (bottom right). The angular power spectra, Cj's, of the CMB-temperature 
anisotropy autocorrelation in different configurations are displayed in figure 4. They have 
been calculated using the Healpix package [64]. The power spectra displayed are defined 
in general for any quantity X as, 

afm = r W r d0 X{9, ^~ X Y lm (0, <t>) (3-1) 
Jo Jo A 

^ m=l 

Q XX = 2l + \ \ a ^m\ > (3-2) 

m=— 1 

with Y[ m being the spherical harmonics. 



3.2 Numerical limitations 

Before discussing the power spectra, a note about numerical limitations needs to be made. 
We can only consider the power spectra up to the multipoles as displayed, because for even 
higher multipoles a much larger resolution in the maps would be needed. Even higher 
resolution would be computationally expensive, up to the point where numerical errors in 
the integration will dominate the anisotropy at such small scales. 

All spectra show an ankle at about the same multipole I ~ 100. The spectrum beyond 
this ankle cannot be trusted, as the ankle is already an aliasing artifact of the resolution 
of choice. When the resolution is increased, the ankle goes down and moves to higher I, 
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for all five configurations displayed. The resolution chosen, 3 x 10 pixels per map, is at 
the limit of the computational resources at hand, with our code as is 3 . 

3.3 Discussion 

By eye one already sees, in figure 3, that the maximum temperature anisotropy is much 
larger than the observed average CMB-temperature anisotropy, notably for larger holes. 
The center of each map points in the special direction of exactly aligned holes. The cold 
spots are directions in which the photons travel through an optimal (read largest) number 
of holes. These are in fact cold, and not hot, spots because travelling through many holes, 
implies crossing many shells of in falling matter as well. What is observed here is actually 
the Rees-Sciama effect. The hot spots correspond to directions in which the photon at t = i 
was in a region that became a hole. I.e., the photon has not travelled through an integer 
number of holes. When crossing an entire hole, a photon passes both regions in which 
the matter, falling outwards to form a shell around the hole, falls against the direction of 
photon and regions in which the matter falls in the same direction as the photon. Any 
Rees-Sciama effect mostly cancels in that case. When the photon is in the middle of a 
region that will become a hole when the matter starts falling, the photon automatically 
will only experience the metric of matter that falls in the direction of the photon path. In 
that hole, the cancellation will not happen, hence the photon relatively gains energy. One 
might wonder wether these hotspots are natural, and if it would not make more sense to 
constrain the number of holes that a photon crosses to be integer, for —0.8 < t < 0. Firstly, 
if the Swiss Cheese is a toy model for structure formation, then there is nothing unnatural 
about the period when the matter perturbations undergo the transition from linear to 
non-linear perturbations. In the linear regime, the potentials do not change. But still, if 
the photon happens to cross such a potential well at the time of transition, the photon 
will become part of a similar hot spot as in this toy model. Secondly, the temperature 
deviation of the hot spots from the average temperature is about the same as the deviation 
of the cold spots. It is likely that removing the hot spots would only change the angular 
correlation by a factor of two. This should be investigated though, keeping in mind that 
by eye one cannot do much statistics on a picture. In this work we choose to stick to one 
well-defined model. 

In each map one also observes the same circular ring of average temperature, centered 
on the map. In the projection and orientation chosen, this ring corresponds to the plane 
surrounding the observer in which the photons always travelled through cheese only. This 
is a direct consequence of the regular distribution of holes chosen. Other observed lines, 
especially visible in the map of smaller holes, are similarly artifacts of special directions in 
the regular distribution of holes. 

Let us consider the angular power spectra, C/'s, of the autocorrelation for these maps, 
in figure 4. The spectra are binned in bins of five multipoles. The Swiss-Cheese spectra 
must be considered as additional secondary anisotropics in the CMB, since the primary 

3 The calculations were performed on an MPI-grid. For the largest holes, the computation of one map 
takes about one hour on 16 cpu's. For the tiniest holes, the computation of a map takes about 30 hours on 
128 cpu's. 
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Figure 5: Left: The maximum luminosity distance (top) and the maximum angular diameter 
distance (bottom) achievable as a function of redshift z, for Swiss-Cheese universes with five different 
hole sizes. These distances would be observed when looking through a series of perfectly aligned 
holes. Note that the maximum lensing effect only marginally changes with hole size. This indicates 
that the ratio between time spent in 'cheese' and time spent in 'holes' along the history of a photon 
path is the quantity that defines the observed distances, as opposed to the physical size of the holes. 
Right: The angular diameter distance as a function of redshift z, comparing the Swiss-Cheese model 
with three different spatially flat FLRW-solutions. 



power spectrum has been ignored. For comparison the power spectrum of the CMB, as 
it would be observed today in the standard ACDM universe, is displayed. For holes of 
radius rhi = 35 Mpc, the secondary anisotropies are of the same magnitude as the CMB 
observed today. This means that a Swiss-Cheese universe with all holes smaller than 
35 Mpc potentially leaves the CMB as observed unaltered. For example, the secondary 
anisotropies for holes of radius 3.5 Mpc are up to two orders of magnitude smaller than the 
observed anisotropies. A decrease of the Rees-Sciama effect with decreasing hole-size was 
also foreseen in Refs. [48, 36]. In a more realistic distribution of holes, the size of holes is 
not a fixed number. In that case, probably a few holes larger than rhi ~ 35 Mpc amongst 
smaller holes, would leave an imprint on the CMB small enough to agree with observations. 

4. Angular diameter distance - redshift relation 

4.1 The same maximum distance for all hole sizes up to 1.75 Gpc 

In the right of figure 5, we repeat for illustration a comparison of cIa(z) for different FLRW- 
cosmologies and one particular Swiss-Cheese cosmology, similar to Ref. [48] . Here we focus 
on the maximum deviation in d& with respect to an FLRW-cosmology with O ma tter = 1. At 
redshift z ~ 1.92, in this Swiss-Cheese cosmology dx corresponds to that in a flat FLRW- 
universe with = 0.4. In Ref. [49] it was argued, however, that a good fit over the whole 
curve to a cosmology with = 0.75, is actually achieved by a Swiss-Cheese universe with 
five holes of radius 250 Mpc between the observer and t = i. Since we are dealing with a toy 
model that probes the ability of structure formation to explain the observed acceleration, 
we prefer to probe the maximum achievable effect at high redshift, over finding the best 
out of a set of perhaps mediocre fits. 
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Figure 6: The full-sky maps of the angular diameter distance at z = 1.92 induced by the Swiss- 
Cheese structure. The quantity shown is dA ^ A ■ Again these are maps in Swiss-Cheese universes 
with holes of size ru = 3.5 Mpc (upper left figure), rhi = 350 Mpc (upper right figure), = 875 
Mpc (lower left figure) and r\ x \ = 1.75 Gpc (lower right figure). In contrast with the redshift maps, 
in the angular diameter distance, the maximum anisotropy does not depend on the size of the holes. 

In the left of figure 5, the maximum g?a is compared for all models. The maximum c?a 
for all models lies in the direction of perfectly aligned holes. In nature such a direction is 
unlikely to exist, but a more natural distribution of holes goes beyond our scope here. The 
maximum effect on <Ia is only marginally dependent on the size of holes. This conclusion 
seems in contradiction with conclusions drawn in Refs. [36, 53], which we address at the 
end of this section. 

4.2 Distribution in different directions 

In figure 6 we show a full sky map of the angular diameter distance at a fixed redshift, 
namely the maximum redshift found in the temperature map of each cosmology at t = i , 
which for each cosmology is close to z = 1.92. This means that in most directions the 
redshift of a photon was slightly less at time t = i, hence the integration was continued 
through cheese only up to the right redshift. Interestingly, the maximum anisotropy in 
angular diameter distance is the same for all models, which was illustrated in the left of 
figure 5. 

In Ref. [51] it was already shown that, for holes of radius 350 Mpc, the average d& 
corresponds to that of a cheese-only universe. In figure 7 we show the angular power 
spectrum of the autocorrelation in the angular diameter map, Cf AdA . Thisfi eure illustrates 
that, even though the maximum effect of the Swiss Cheese on <Ja is about the same for 
all hole sizes, the number of special directions in which there is a high effect (or anti- 
effect), decreases with decreasing hole size. The smaller the size of holes is, the smaller the 
standard deviations in d,A will be. This is also apparent in figure 8, where the distributions 
of redshift z and the angular diameter dx are compared for Swiss-Cheese universes with two 
different hole sizes. The standard deviation in redshift goes down from = 2.6 x 10 -4 
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Figure 7: The C? ac1a, s for different cosmolo gics, in bins of 5 multipoles. Even though the maximum 
observed d A is independent of the size of holes, the average anisotropy goes down with decreasing 
hole size. This is because the number of special directions decreases with decreasing hole size. 




0.3415 



0.341 




Figure 8: Left: The distribution of redshift of photons in all directions, at fixed time. Clearly, 
for smaller holes, the variance in redshift decreases. Right: The distribution of angular diame- 
ter distance of photons in all directions, at fixed redshift. Again, for smaller holes, the variance 
decreases. 

for rhi = 350 Mpc to a_^_ = 3.8 x 1CP 5 for r^i = 35 Mpc. Similarly for the angular diameter 
distance, it goes down from a^ A = 1.3 x 10~ 2 to Od A = 4.3 x 10 -3 . 

4.3 Dependence on the size of holes 

In Refs. [36, 53] it was found that the effect of one hole on cZa, with respect to a cheese-only 
passage, goes as Ag?a = ^a° 1c ~~ d c ^ eese ~ r hi/-^H' w ith -^h the Hubble radius. Increasing 
the number of holes with N^i ~ Rn/r^i, this would mean that the maximum effect scales as 
^A,max ~ r hi/-^H- As total lensing effect is crucial to the success of a Swiss-Cheese model, 
this contradictions deserves more attention. The difference in size-dependence may be an 
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Figure 9: An ambiguity when speaking of the effect of the size of holes on d A - Left: Two sketches 
of methods comparing hole sizes. The first method considers only distance traveled through a hole, 
and compares c?a at the exit of that hole to that of an cheese-only photon. The second method 
considers one central void embedded in cheese, and compares cIa at the end of the same total 
comoving distance through a hole and cheese. Right: Using the first comparison method we recover 
the behaviour that is found in Refs. [36, 53], Ad A oc rjj (red, solid). Using the second method, we 
find a different behaviour which is closer to Ad A oc r^i (blue, dashed). 



artefact of the modeling. However we find a difference at the analytical level, before even 
deciding which model to take. Let us explain why we come to a different conclusion. 

The calculations in Refs. [36, 53] consider one hole at a time, and one hole only. The 
effect of several holes does however not grow linearly with the number of holes. It may 
happen to do so for a particular modeling, but in general it will not. In figure 9 we show 
the results of two different approaches. The first approach considers only one hole. If one 
starts an integration when entering a hole, and stops when exiting the hole in order to 
compare d A with the value found when travelling the same distance through cheese, we 
recover in fact the foreseen effect of Ad\ oc rjj. In the second approach, we consider a 
photon passing one hole along a fixed distance, where all space outside the single hole is 
necessarily cheese. In this case we do not find the same behaviour, but the size dependende 
seems to be more like Ad a oc r^\- 

The reason for this discrepancy is that, after passing a void, the quantity (see 
equations (2.20,2.22)) keeps a deviation from the same quantity for a cheese-only photon, as 
shown in figure 10. Hence, if one were to consider the part of the integration, starting only 
beyond the patch, the initial conditions for the cheese-only photon and the photon that 
passed through the hole are different. The quantity can be interpreted as a measure 
of divergence / convergence of the beam. Changing this quantity necessarily changes the 
final beamsize, even when no more holes are on the way. There is no reason to believe 
that this quantity obtains no change from a hole, as 'the magnification of the beam, is 
proportional to the integrated column density along the beam's path' [51]. The integration 
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Figure 10: Left: A comparison of the dd^/dr along the trajectories of two photons. Time increases 
from right to left, and the initial conditions are chosen such that an observer today (at r — —0.21) 
sees a source through cheese only (red solid line) and a source through one central hole, with 
rhi = 350 Mpc, centered between the observer and t — t (green dashed line). Both sources are 
observed at the same solid angle, which translates to d&^/dr at the observer. They have however 
a different size, which translates to d&^/dr at the source. That means that the source at r = 0.21 
behind the hole must have a larger area in order to be observed under the same solid angle, and 
hence has a larger angular diameter distance. The key point in this figure is the difference in 
ddx/dr at both sides of the hole. Right: A[ddA/dr) shows the suppressing effect the hole has on 
the derivative of the area of the beam, causing the source to be observed under smaller solid angle 
than a source of the same size seen through cheese only. Hence, once a photon passes a hole, the 
lensing effect of the hole continues playing a role beyond the hole, by having changed the initial 
conditions of the integration beyond the hole. 



is performed over the past light cone of a photon. The photon spends more time in the 
void of a hole than it does in the shell of a hole, therefore there is no cancellation of the 
defocussing effect when leaving a hole. In fact, a closer look at Eqs. (2.7-2.10,2.21,2.22) for 
the special case of a radial photon, such that c<h = £ = 0, shows us, 



d 2 VA 
dr 2 



2 r- Y'(r,t) 2 p(r,t) 

3 W{r) 2 



+ 



dr 



Y"(r,t) _W\r) Y'(r,t) 
Y'(r, t) ~ W{r) + W(r) 



VAF 1 (r,t) + d ^F 2 (r,t). 
dr 



(4.1) 



Here we used, 



d?_ _ UX\ 2 d 2 + £\d 
dr 2 \dr J d\ 2 dr 2 d\ ' 



(4.2) 



in order to be able to unambiguously compare photons along different geodesies, with 
different afine parameters A but with the same coordinates r and t. Now let us compare 
two photons, that entered a patch with the same initial conditions, and let one of them see 
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only cheese. At the exit of the hole, the difference in rf ^2 is, 



A 



d 2 VA 
dr 2 



A 



F 1 (r hh t) + A 



dy/A 
dr 



F 2 (r hh t) 



(4.3) 



where in we use the notation A[Q] = Q hoie photon (nil) - Qcheese photon (ill)- 

First consider Fi(rh\,t) and i*2( r hi> *)• These are functions of background quantities 
only, evaluated at r = r^i where all background quantities are supposed to continuously 
match from the hole to the cheese. Hence, these functions are the same for all photons, 
no matter whether they saw a hole or cheese only, no matter what initial conditions their 



beam had. By construction we have A 



/ 0, since that is the whole purpose of 



the Swiss-Cheese structure. Therefore, beyond the hole, one expects that at least either 



A 
ho 



d-jA 



dr 



/ or A 



d 2 y/A 



W 1 
the hole. 



/ if A 



7^ 0, proving that the defocussing effect of the 

e carries beyond 

The reader should keep in mind that the defocussing effect is of importance when one 
considers the maximum effect one hole can have on d^, even beyond the passage of that 
particular hole. Remember that we are considering the divergence of a beam, after having 
passed a hole, which is the derivative of the angular diameter distance. However, as it 
has been shown in Refs. [51] that the average effect of one hole on the angular diameter 
distance disappears when averaged over angles of incidence, it is likely that the average 
effect of one hole on the divergence of the beam (which, as we just discussed, influences 
the angular diameter distance beyond the hole) will also disappear on average. 

A final remark to make here, is a reminder that we are dealing with a toy model. The 
model is constructed manually and ideally. The defocussing by the hole may be an artefact 
of the modeling, and may disappear in a more realistic model. We showed that in idealized 
models it is most likely that the divergence of a beam is altered by a hole, but this gives 
no guarantee for more realistic models. 



5. Discussion and conclusion 



For the first time we have performed a full sky simulation of the CMB in Swiss-Cheese 
universes, in realizations with different hole sizes. We have shown that if all holes have 
a radius larger than r^i = 35 Mpc, the Swiss-Cheese model can be ruled out on basis of 
the observed CMB. One could try to change the density profile of the holes, but in order 
to leave the CMB intact, that new profile would necessarily be closer to FLRW, hence no 
longer reproducing the wanted angular diameter distance effect. 

Another option to try to save the model is to decrease the size of holes. We find that 
for small holes, the maximum angular diameter effect remains the same as for larger holes. 
However, the probability that special directions exist, in which one sees through the centers 
of a series of perfectly aligned holes, decreases with decreasing hole size, at least in the 
case of holes on a lattice, as in this work. 

In Refs. [44, 51] it was argued that (spherical) mass compensated holes in a Swiss- 
Cheese universe cannot make the average angular diameter distance depart from the EdS- 
case, when the observer is outside the hole. Not surprisingly, we confirm this result. 
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Moreover, we have shown that even though the average d\ does not change, the anisotropy 
in CMB-temperature due to the same hole does in fact differ significantly from the EdS- 
case, which in principle should not depend on the sphericalness of a hole, but on the mere 
fact that the observer sees both through and past the side of the hole, inevitably leaving a 
Rees-Sciama imprint on the CMB. 

These findings together lead to the conclusion that the Swiss-Cheese universe, as con- 
sidered here, with mass compensated holes, all of them larger than = 35 Mpc, spherical 
or not, is ruled out. Based on the findings in Ref. [51], and the findings in this work with 
respect to the variance in d\ in different directions that decreases with decreasing hole size, 
also Swiss-Cheese universes with spherical mass compensated holes of radius smaller than 
rhi = 35 Mpc are ruled out. 

This conclusion applies to all Swiss-Cheese-like models in which the universe is satu- 
rated with mass compensated voids, and that try to explain the angular diameter distance 
at high redshift. We did not consider the case where our local Hubble volume contains 
just a small number of voids, without being saturated with them. Such a scenario remains 
interesting in order to explain anomalies in the CMB, such as the cold spot [61]. How- 
ever, such a scenario is in principle unrelated to possible alternative explanations for the 
observed accelerated expansion. 

We have shown that the maximum effect of holes on g?a has weaker, if not zero, 
dependence on the size of the holes than previously claimed in the literature, due to the 
defocussing effect a hole has, which carries influence beyond the hole. We assessed the 
importance of taking into account this defocussing, when addressing the effect of different 
hole sizes. This effect could be due to the modeling of the density profile. The model 
dependence should be more carefully assessed, before definitive conclusions can be drawn. 
It remains to be investigated as well, what the effect would be of non-spherical holes, with 
a typical size < 35 Mpc, densely spread in such a manner that locally the universe is 
nowhere described by the cheese metric, but globally remains FLRW. Another addition to 
subject to research, is to see the effect of virialization of the mass shells, which, as pointed 
out in Ref. [36], may increase the effect on d\, and at the same time should decrease the 
Rees-Sciama effect, hence decreasing the potential trouble with the CMB. 
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